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ABSTRACT:  Gel  systems  based  on  self-assembled,  amphiphilic 
ABA  triblock  copolymers  in  midblock-selective  solvent  form  sta¬ 
ble,  spatially  extended  networks  with  controllable  morphology 
and  tunable  viscoelastic  behavior.  In  this  work,  we  systemati¬ 
cally  evaluate  the  mechanical  properties  of  these  gels  using  mor¬ 
phology  calculations,  and  a  nonequilibrium  oscillatory  shear 
technique  based  on  the  dissipative  particle  dynamics  (DPD) 
method.  Our  simulations  demonstrate  that  low  molecular 
weight  triblock  copolymers  with  incompatible  blocks  self-assem- 
ble  into  micelles  connected  with  bridges  and  loop-like  chains 
comprised  of  the  solvent-selective  polymer  midblocks.  The  frac¬ 
tion  of  bridges,  </>b,  generally  increases  with  increasing  relative 
volume  of  the  midblock,  x,  defined  as  the  ratio  of  midblock  and 
endblock  volumes  (x  =  ^).  For  our  model,  (£b  reaches  a  plateau 
at  approximately  x  >  9  for  a  strongly  selective  solvent.  At  this 
limit,  the  value  of  4>b  increases  from  0.40  to  about  0.66  as  the 


copolymer  concentration,  c,  increases  from  0.2  to  0.5;  however, 
this  increase  is  less  significant  at  higher  concentrations.  The 
elastic  response  of  the  gel  studied  here  is  comparable  with  the 
Rouse  modulus.  The  elastic  modulus  increases  with  polymer 
concentration,  and  it  exhibits  a  broad  peak  within  6  <  x  <  12. 
Finally,  we  present  an  approximate  method  to  predict  the  elas¬ 
tic  modulus  of  unentangled  ABA  triblock  copolymers  based 
solely  on  the  morphology  of  the  micellar  gel,  which  can  be 
gleaned  from  equilibrium  DPD  simulations.  We  demonstrate 
that  our  simulation  results  are  in  good  qualitative  agreement 
with  other  theoretical  predictions  and  experimental  data.  ©  2009 
Wiley  Periodicals,  Inc.*  J  Polym  Sci  Part  B:  Polym  Phys  48:  15- 
25,  2010 

KEYWORDS:  block  copolymers;  DPD  simulation;  gels;  morphol¬ 
ogy;  viscoelastic  properties 


INTRODUCTION  Polymer  gels  based  on  blends  of  synthetic 
block  copolymers  in  solvent  exhibit  thermal  and  environmental 
stability  and  possess  controllable  morphological  and  mechani¬ 
cal  properties.  Amphiphilic  ABA  triblock  copolymers  are  one 
example  of  a  polymer  gel  that  has  been  widely  used  in  many 
practical  applications.  ABA  triblock  copolymers  undergo  ther- 
moreversible  microphase  separation  in  midblock-selective  sol¬ 
vent  and  form  spatially  extended  networks.1'2  The  microstruc¬ 
ture  of  these  gels  depends  on  several  factors  such  as 
temperature,  copolymer  concentration,  block  architecture,  and 
relative  size  of  the  hydrophobic  and  hydrophilic  blocks.  The 
polymer  networks  that  form  consist  of  physical  crosslinks  com¬ 
prising  of  solvent-incompatible  endblocks  (micelles]  connected 
by  the  midblocks  (bridges].  The  micelles  are  surrounded  by  a 
corona  consisting  of  midblocks  (loops],  formed  as  a  result  of 
the  close  association  of  the  endblocks.  Some  endblocks  may 
not  assemble  into  micelles,  producing  so-called  tails  and  free- 
chains,  where  a  tail  has  only  one  endblock  associated  with  the 
network,  whereas  neither  endblocks  of  a  free-chain  are  part  of 


the  network.  A  schematic  of  such  a  polymer  network  is  given 
in  Figure  1.  For  unentangled  triblock  copolymers,  the  bridges 
can  be  considered  as  active  chains  of  the  network  and  are 
treated  as  elastic  springs,  whereas  contributions  to  the  gel  net¬ 
work  elasticity  from  loops,  tails,  and  free-chains  can 
be  considered  negligible.3'4  For  that  reason,  the  bridge-to-loop 
ratio  or  bridge  fraction  is  one  of  the  key  quantities  that  govern 
the  mechanical  properties  of  a  gel.  The  bridge  fraction  depends 
on  several  factors,  such  as  copolymer  concentration,  degree  of 
block  incompatibility,  and  the  relative  size  of  the  blocks.  Exper¬ 
imental  observations  have  revealed  that  the  bridge  fraction 
increases  with  relative  size  of  the  midblock  (with  respect  to 
the  endblocks]  and  copolymer  concentration.5-7 

The  network  resulting  from  bridging  strongly  influences  sev¬ 
eral  mechanical  properties  that  are  often  used  to  character¬ 
ize  the  gels.  For  example,  when  a  linear  viscoelastic  material 
is  subjected  to  oscillatory  strains,  y,  of  frequency,  co,  the 
stress  response,  cr,  is  necessarily  cyclic  and  can  be  written 
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FIGURE  1  Schematic  representation  of  a  polymer  network 
formed  from  an  ABA  triblock  copolymer  in  a  midblock-selective 
solvent. 


as:  o[co)  =  (G/(co)sin(cot]+  G"(<x>]cos(cut]],  where  Am  is 
the  amplitude.  The  storage  modulus,  G'[co ],  characterizes  the 
material's  ability  to  store  energy  and  the  loss  modulus, 
G"[co),  characterizes  the  amount  of  energy  lost  through  vis¬ 
cous  processes.  For  gels,  G'[co)  often  is  only  weakly  depend¬ 
ent  on  co  over  a  broad  range  and  is  several  orders  of  magni¬ 
tude  larger  than  G"[co).  Thus,  a  reasonable  assumption  is  to 
define  the  elastic  response  of  a  gel  by  referring  to  the  pla¬ 
teau  value  of  G1  or  the  equilibrium  elastic  modulus,  Ge8  Ge 
for  an  ABA  gel  made  from  unentangled  chains  in  the  semidi¬ 
lute  regime  in  a  selective  solvent  is  comparable  with  modu¬ 
lus  of  Gaussian  network  in  a  good  solvent  Ge>8,9  expressed  as 


^Rouse 


(i) 


where  C  is  a  mass/volume  concentration,  Mw  is  the  molar 
mass  of  the  triblock  copolymer,  T  is  the  temperature,  R  is 
the  universal  gas  constant,  and  0b  is  the  previously  defined 
bridge  fraction.  For  copolymers  with  higher  molecular 
weight,  entanglements  would  give  rise  to  a  large  additional 
contribution  to  the  elastic  modulus,  and  the  value  of  Ge 
would  be  an  order  of  magnitude  higher  than  that  predicted 
with  eq  1.  According  to  classical  polymer  theory,  Ge  for 
entangled  polymers  is  proportional  to  1/Nec23,  where  Ne  is 
the  number  of  concentration  blobs  per  entanglement  with 
size,  £  =  c-1^,  c  is  the  volume  fraction  of  polymer,  and  v 
is  the  Flory  exponent  (v  =  0.588]. 10,11 

Several  studies  using  molecular  simulation  have  corrobo¬ 
rated  and  expanded  the  theoretical  picture  of  the  morphol¬ 
ogy  and  mechanical  properties  of  gels  formed  from  ABA  tri¬ 
block  copolymers  in  midblock-selective  solvent.12-16  Khalatur 
et  al. 17,18  studied  the  micellar  organization  and  rheology  of 
the  triblock  gels  and  found  that  equilibrium  microstructures 
of  these  gels  depends  on  the  chain  rigidity,  and  that  their 
rheological  properties  are  governed  by  the  spatial  distribu¬ 
tion  of  bridges  between  micelles. 

Moreover,  there  have  been  several  theoretical  and  computa¬ 
tional  efforts  undertaken  to  gain  a  physical  understanding  of 
the  role  of  the  bridge  fraction  on  the  network.  By  using  self- 
consistent  field  theory,  Zhulina  and  Halperin19  have  demon¬ 


strated  that  the  bridge  fraction  is  proportional  to  the  degree 
of  incompatibility,  </>b  oc  (xM)-1^9.  Matsen  and  Schick20  have 
shown  that  </>b  decreases  slowly  with  increasing  polymer 
length,  M,  increasing  Flory-Huggins  parameter,  y,  and 
decreasing  endblock  fraction  in  the  strongly  segregated  la¬ 
mellar  phase.  Although  Matsen  and  Schick  found  that  the 
bridge  fraction  was  typically  equal  to  0.40,  Jones  et  al.21 
found  the  bridge  fraction  to  be  approximately  equal  to  0.63 
for  the  strong  segregation  limit.  An  increase  of  cj)b  with  co¬ 
polymer  concentration  has  also  been  reported  by  Monte 
Carlo  simulation  studies,4,22-24  where  Karatasos  et  al.23  esti¬ 
mated  the  bridge  fraction  to  be  in  the  range  of  0.37-0.50, 
with  increasing  molecular  weight  of  the  triblock.  Nguyen- 
Misra  and  Mattice  have  demonstrated  that  the  bridge  frac¬ 
tion  increases  with  copolymer  concentration  and  block 
incompatibility.  The  same  authors  observed  that  the  bridge 
fraction  was  slightly  higher  with  an  increase  of  the  relative 
midblock  size,  but  found  this  increase  to  be  insignificant.24 
These  theoretical  findings  have  been  partially  confirmed  by 
experimental  work.7,25  Watanabe  et  al.7  demonstrated  that 
the  bridge  fraction  increases  from  0.2  to  0.4  with  increasing 
copolymer  concentration  from  20  to  50  weight  percent. 
Takano  et  al.  obtained  a  bridge  fraction  equal  to  0.93,  which 
is  higher  than  theoretical  predictions  and  other  experimental 
work.  They  attributed  this  finding  to  the  possible  existence 
of  interlocked  entanglements.26 

Despite  the  significant  impact  that  these  computational  and 
theoretical  studies  have  made  to  our  understanding  of  the 
structure  and  mechanical  properties  of  triblock  copolymers, 
there  are  some  considerable  shortcomings.  Micellar  phases 
are  characterized  by  fluctuation  effects;  hence,  they  may  not 
lend  themselves  to  an  accurate  implementation  of  methods 
such  as  self-consistent  field  theories.  Furthermore,  the  Monte 
Carlo  approach  does  not  lend  itself  to  the  simulation  of  the 
dynamics  of  fluids.  Moreover,  the  impact  of  chemical  and 
compositional  variables  on  the  bridge-to-loop  ratio  of  tri¬ 
block  gels  has  yet  to  be  systematically  studied  by  a  particle- 
based  dynamics  approach  at  either  the  molecular  or  meso- 
scale.  To  address  this  shortcoming,  we  present  a  dissipative 
particle  dynamics  (DPD]  simulation  study  of  the  ABA  tri¬ 
block  copolymer  in  a  B-selective  solvent.  In  this  study,  we 
attempt  to  create,  through  computational  modeling,  a  predic¬ 
tive  tool  that  answers  questions  regarding  the  impact  of  con¬ 
centration  and  molecular  architecture  on  the  morphology,  in 
particular  the  bridge  fraction,  and  viscoelastic  properties  of 
ABA  gels.  Understanding  such  structure-property  relation¬ 
ships  will  enable  the  development  of  thermoreversible  gels 
with  the  desired  morphological  and  mechanical  properties 
that  meet  the  specific  application  needs. 

MODEL  AND  SIMULATION  METHODS 

Generally  speaking,  atomistic  molecular  dynamics  is  suitable 
up  to  10  nm  and  several  nanoseconds,  making  it  difficult  to 
perform  simulations  of  large  polymeric  systems  on  compara¬ 
ble  timescales.  In  this  study,  we  have  exploited  the  DPD  sim¬ 
ulation  technique,  consisting  of  coarse-grained  particles  that 
represent  clusters  of  molecules  rather  than  individual 
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atoms.27'28  The  DPD  method  is  a  mesoscale  simulation  tech¬ 
nique  that  operates  at  time  and  length  scales  larger  than 
those  of  traditional  molecular  dynamics,  whereas  providing 
the  ability  to  probe  phenomena  that  is  inaccessible  to  contin¬ 
uum  modeling. 

Dissipative  Particle  Dynamics 

A  DPD  system  is  composed  of  soft  particles,  each  represent¬ 
ing  a  region  of  fluid,  moving  according  to  Newton's  equa¬ 
tions  of  motion  continuously  in  space  and  discreetly  in  time. 
In  a  DPD  simulation  of  polymers,  the  chain  is  modeled  as  a 
collection  of  point  particles  that  represent  lumps  of  the  chain 
containing  several  segments.  DPD  particles  are  defined  by 
mass,  rrij,  position,  r/;  and  velocity,  vz,  and  interact  with  each 
other  via  a  pairwise,  two-body,  short-ranged  force,  F,  that  is 
written  as  the  sum  of  a  conservative  force,  Fc,  dissipative 
force,  Fd,  and  random  force,  FR,  as  follows: 

F‘  =  EFEEF°+EF*  pj 

i+i  i+i 

where  Fc  includes  a  soft  repulsion  force,  FCr,  acting  between 
two  particles  and  a  Finite  Extensible  Nonlinear  Elastic  (FENE] 
force,  Ffene,  acting  between  adjacent  particles  in  a  polymer 
chain.  The  DPD  polymer  chains  defined  in  this  manner  are 
flexible,  because  no  additional  constraints,  such  as  bond  bend¬ 
ing  or  bond  torsion,  are  included.  FCr  and  feene  are  given  by 


In  eqs  3  and  4,  azy  is  the  maximum  repulsion  between  parti¬ 
cle  /  and  particle  j,  r,y  =  |rz-  —  r;  |,  ezy  =  (rz  —  r7)  /rzy  is  the  unit 
vector  from  the  jth  particle  to  the  ith  particle,  rc  is  the  cutoff 
radius,  H  is  the  spring  constant,  /  represents  the  maximum 
extension  length  of  a  polymer  bond,  and  r  is  the  vector  that 
connects  pairs  of  adjacent  particles  in  a  polymer  chain.  The 
remaining  two  forces,  F°  and  FR,  are  given  by 

F,z  =  -ycoD  (rv)  (e,j-Vij)  ev,  (5) 


and 


F|*  =  (7ft)R 


VL< 
Va t 


(6) 


where  coD(r)  and  coR(r)  are  weight  functions,  y  is  the  friction 
coefficient,  a  is  the  noise  amplitude,  vzy  =  vz  —  Vy,  ^  is  the 
Gaussian  random  number  with  zero  mean  and  unit  variance 
that  is  chosen  independently  for  each  pair  of  interacting 
particles,  and  At  is  the  time  step. 

Espanol  and  Warren29  showed  that  the  system  samples 
the  canonical  ensemble  and  obeys  the  fluctuation-dissipation 
theorem  (in  the  limit  of  At^O]  if  the  following  relations 
hold. 


coD(r)  =  [coR(r)]2, 


(7) 


and 


<72  =  2  ykBT,  (8] 

where  T  is  the  temperature,  and  kB  is  Boltzmann's  constant. 
coR(r]  is  typically  chosen29  as 


coR(r) 


1  -  £  for  r  <  rc 
0  for  r  >  rc  ' 


(9) 


For  a  more  detailed  description  of  the  DPD  method,  see  the 
original  articles  by  Hoogerbrugge  and  Koelman.27'28 

Coarse-Grained  Model  of  A1BXA1  Triblock  Copolymer 

The  three-dimensional  system  consists  of  a  total  of  Nt  DPD 
particles  of  triblock  copolymers  connected  by  a  FENE  force, 
plus  solvent  in  a  cubic  domain  with  periodic  boundaries  in 
all  directions.  The  system  has  N  chains  of  an  A1BXA1  triblock 
copolymer  comprised  of  particles  of  types  "A"  and  "B"  of 
total  length  M  immersed  in  a  midblock-selective  solvent 
whose  particles  are  of  type  "S."  Each  DPD  particle  represents 
roughly  the  same  volume.  Each  chain  has  two  endblocks  rep¬ 
resented  by  a  single  particle  of  type  'A"  and  a  midblock  com¬ 
prised  of  x  particles  of  type  “B."  See  Figure  2  for  a  schematic 
of  an  A1B3A1  triblock  copolymer. 

The  size  of  the  DPD  particle  in  the  simulated  chain  should 
be  larger  than  or  equal  to  a  Kuhn's  segment  of  the  polymer 
to  ensure  the  flexibility  of  the  polymer  on  the  mesoscopic 
scale.30 

The  number  of  triblock  chains  is  chosen  from 


(10) 


where  M  =  x  +  2  and  c  is  the  volume  fraction  of  the  triblock 
in  the  solution,  and  x  is  the  number  of  particles  making  up 
the  midblock.  The  endblocks  for  all  triblocks  are  modeled  as 
a  single  particle  to  avoid  any  occurrence  of  the  chain  collaps¬ 
ing  on  itself  before  any  association  with  other  endblocks 
takes  place. 

The  physical  size  of  the  interaction  radius  for  the  DPD  parti¬ 
cle,  rc,  depends  on  the  molar  volume  of  the  endblock31  and 
may  be  estimated  from 


_  _  s/Pdpdva 

c  V  Na  ' 


(ii) 


where  pDPD  =  ^  is  the  DPD  particle  density,  Vt  is  the  total 
volume  of  the  system,  vA  is  the  molar  volume  of  the  end- 
block,  and  Na  is  Avogadro's  number.  The  number  of  mid¬ 
block  particles,  x,  is  evaluated  from  the  volume  fraction  of 
the  endblock  in  the  dry  copolymer,  fA>  written  as 


2(1 -A) 

/a  ' 


(12) 
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FIGURE  2  Schematic  representation  of  an  A-iBsA-i  triblock  co¬ 
polymer  mapped  onto  a  DPD  model.  The  poly(styrene- 
isoprene-styrene)  copolymer  is  chosen  as  an  example. 


delineating  the  regimes  and  frequencies  at  which  various 
effects  (e.g.,  polymer  chain  length}  are  manifested.  Because 
of  the  soft  nature  of  the  DPD  potentials  used  in  this  study 
polymer  chains  can  pass  freely  through  each  other,  thus  our 
simulations  cannot  mimic  chain  uncrossability.32  The  visco¬ 
elastic  properties,  such  as  storage,  G'[co ),  and  loss  moduli, 
G"[(d),  can  be  calculated  using  a  nonequilibrium  oscillatory 
shear  technique,33'34  entailing  a  simulation  with  an  addi¬ 
tional  force  in  the  shear  direction,  along  with  time-depend¬ 
ent  Lees-Edwards  boundary  conditions.35  For  oscillatory 
shear  imposed  in  the  x-y  plane,  the  equation  of  motion  for 
the  particle  velocities  becomes34 


dvix  d2y(t) 


(13) 


The  oscillatory  strain  has  been  taken  as  y(t]  =  Am  sin  [cot], 
where  and  co  are  values  chosen  for  the  amplitude  and 
frequency,  respectively.  We  have  used  =  0.1  that  corre¬ 
sponds  to  the  linear  rheology  regime  and  an  oscillation  pe¬ 
riod,  p  =  2n/co  from  101  to  104. 

The  equilibrium  DPD  approach  allows  us  to  define  the  vari¬ 
ous  conformations  of  the  midblocks  (bridges,  loops,  and  dan¬ 
gling  tails]  and  precisely  calculate  their  quantities.  After¬ 
wards,  we  can  correlate  the  structural  properties  of  the 
triblock  gel  with  the  elastic  modulus  obtained  from  the  oscil¬ 
latory  shear  simulations. 


Note  that  because  the  endblock  of  any  length  is  represented 
by  one  DPD  particle,  x  corresponds  to  the  relative  volume  of 
the  midblock  in  the  triblock  copolymer:  x  =  VB/VA,  where  VA 
and  VB  are  volumes  of  the  midblock  and  endblock, 
respectively. 


Assuming  that  the  crosslinks  [micelles]  are  fixed  and  moving 
in  an  affine  manner,  the  equilibrium  elastic  modulus,  Ge,  can 
be  evaluatedas  a  sum  of  the  contributions  of  the  different 
conformations  of  the  midblock,7  as  follows: 

Ge  =  Vb^b  +  v\G\  +  vtGt  [14] 


Viscoelastic  Property  Calculations 

Simulations  have  been  used  to  determine  the  parameters 
that  govern  the  viscoelastic  behavior  of  polymer  systems, 


where  p  is  the  copolymer  density,  vb,  Vi,  and  vt  are  the  num¬ 
ber  densities  of  bridges,  loops,  and  dangling  tails,  respec¬ 
tively,  and  where  Gh,  G\,  and  Gt  are  the  elastic  modulus  of  a 


FIGURE  3  Snapshots  from  a  DPD  simulation  of  our  model  A1B6A1  copolymer  for  the  system  size  of  15000  DPD  particles,  where 
micelles  comprised  of  type  "A"  particles  are  colored  blue,  and  bridge  and  loop  conformations  of  the  midblocks  of  type  "B"  par¬ 
ticles  are  colored  red  and  brown,  respectively.  For  visual  clarity,  solvent  particles  are  not  shown.  The  equilibrated  structures  are 
for  (a)  c  =  0.2  and  (b)  c  =  0.4. 
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TABLE  1  Repulsive  Parameters  Between  a  Pair  of  DPD  Particles 

DPD  Pair  S-S  A-A  B-B  A-B  A-S  B-S 

a  ij  50.0  50.0  50.0  66.0  105.0  50.0 

single  bridge,  loop,  or  dangling  tail,  respectively.  In  the  high- 
segregation  regime,  the  number  of  tails  is  negligible  and  the 
midblocks  adopt  either  a  bridge  or  loop  conformation.  By 
using  eq  10  and  evaluating  the  number  of  bridges,  Nh  = 
iV</>b(x,  c),  and  the  number  of  loops,  iV/  =  N[1  —  cpb[x,  c]),  we 
can  determine  the  equilibrium  elastic  modulus  in  DPD  units 
as: 

Ge  =  Pdpd  [GtAO,  c)  +  G/(  1  -  *,(*,  c))]  (15) 

where  0b(x,  c)  is  the  bridge  fraction  and  is  a  function  of  the 
relative  block  size  and  copolymer  concentration.  With  </>b(x, 
c]  and  Ge  evaluated  from  the  equilibrium  and  oscillatory 
shear  simulations,  respectively,  we  can  estimate  the  contribu¬ 
tions  of  a  single  bridge  or  loop  to  the  elastic  modulus. 

Simulation  Parameters 

DPD-reduced  units  are  adopted  for  the  convenient  expres¬ 
sion  of  parameters  and  are  taken  as  length  in  rc,  energy  in 
kBT/rcl  mass  in  particle  mass  m,  and  time  in 
t  =  rc^/m/kBT.31  Therefore,  the  actual  value  of  the  moduli 
are  calculated  from  the  moduli  expressed  in  DPD  units,  GDPD 
as 

G  =  GDPDkBT/(rc)3  (16) 

We  choose  standard  values  of  m  =  kB T  =  rc  =  1,  time  step 
At  =  0.01  and  overall  particle  density  pDPD  =  3.  H  and  /  in 
eq  4  for  the  triblock  are  1  and  3,  respectively,  which  gives 
the  length  of  a  single  bead-spring  unit  equal  to  one  Kuhn's 
segment,  Z?Dpd  =  =  1-36 

We  have  chosen  the  repulsive  parameter  for  like-particles, 
an,  for  all  types  of  particles  to  be  equal  to  50.  Groot  and 
Warren  determined  an  by  matching  the  dimensionless  com¬ 
pressibility  for  a  real  system,  kT1,  and  the  DPD  fluid  accord¬ 
ing  to  the  scaling  rule  an  =  kB T k2(XNpm~ 1  >  where  Nm  is  the 
level  of  coarse-graining  and  a  is  determined  to  be  0.101  ± 
0.00 1.31,37  For  water  at  pDPD  =  3,  kBT  =  1  and  under  stand¬ 
ard  conditions  (k_1  =  15.98],  Trohimov  estimated  the  maxi¬ 
mum  level  of  coarse-graining  should  be  iVJJJax  =  10  water 
molecules  to  represent  one  DPD  particle.38  Having  Nm  larger 
than  iVJJJax  would  give  an  >  250  and  cause  an  artificial  solidi¬ 
fication  of  the  liquid  system.38'39  Thus  by  modeling  tens  of 
polymer  units  as  a  single  DPD  particle,  we  are  greatly 
exceeding  the  value  of  Conversely,  Nakamura  and 

Tamura40  have  shown  that  the  dependence  of  the  morphol¬ 
ogy  of  an  A-B  type  chain  in  an  A-selective  solvent  on  ah  is 
weaker  than  the  dependence  on  a,y.  It  is  not  always  possible 
to  correctly  map  an  atomistic  simulation  to  a  mesoscale  rep¬ 
resentation;41  therefore,  the  choice  of  az/  =  50  seems  to  be 
reasonable. 


The  repulsive  parameters  for  unlike-particles  a zy  are  deter¬ 
mined  according  to  a  linear  relationship  with  the  Flory-Hug- 
gins  parameter37: 

aij  «  an  +  3.27/;y.  [17] 

where  Xij  parameters  are  mapped  from  long  polymers  to  the 
shorter  DPD  chains  while  preserving  %M  =  const.  We  have 
chosen  the  excess  repulsion  Aa  =  3.27xzy  for  an  A-B  pair 
equal  to  55.  This  value  corresponds  to  the  strong  segregation 
regime,  %M  ^  85  —  340,  where  M  changes  from  5  to  20. 
Interaction  parameters  between  blocks  of  the  copolymer  and 
solvent  are  chosen  such  that  Xas  ~  5  and  xbs  =  0,  respec¬ 
tively,  to  model  a  perfectly  selective  solvent,  i.e.,  nonsolvent 
for  the  A-block  and  good  solvent  for  the  B-block. 

The  maximum  repulsive  force  parameters  used  in  this  study 
are  summarized  in  Table  1. 

Methodology 

This  study  investigates  the  role  of  the  relative  block  size  and 
triblock  copolymer  concentration  on  the  structure  and  struc¬ 
ture-related  mechanical  properties  of  the  polymeric  gels  that 
form.  To  probe  the  effect  of  concentration  and  relative  block 
size,  we  considered  a  range  of  values  for  both  c  (0.2,  0.3,  0.4, 
and  0.5)  and  *.3.6.9,12,15,18 

Equilibrium  configurations  were  generated  by  randomly 
placing  polymer  chains  and  solvent  particles  into  the  simula¬ 
tion  box,  followed  by  slowly  adjusting  the  periodic  box  size 
to  the  desired  pressure.  This  was  followed  by  a  constant 
temperature  simulation  for  1000  steps.  An  equilibrium  state 
was  determined  by  repeating  this  procedure  six  times  until 
the  system  pressure  was  stable  within  1%.  Ten  parallel  sim¬ 
ulations  were  run  for  each  set  of  initial  conditions  to  mini¬ 
mize  statistical  uncertainty.  After  equilibration  of  the  struc¬ 
ture,  a  series  of  stress  trajectories  was  generated  by 
imposing  the  oscillatory  shear  conditions  with  0.1%  strain. 
Because  DPD  simulations  are  generally  able  to  access  only 
frequencies  in  the  megahertz  regime,  dynamic  moduli  were 
linearly  extrapolated.  Note  that  numerical  errors  tend  to  be 
higher  at  low  frequencies  requiring  averaging  over  more 
stress  trajectories  to  minimize  these  errors. 

We  examined  the  microstructure  of  each  equilibrium  state  to 
determine  the  conformation  of  each  midblock  in  the  gel 
using  the  following  scheme.  First,  we  identified  all  crosslinks 
or  micelles.  A  micelle  was  defined  as  two  or  more  endblocks 
in  close  proximity,  such  that  any  two  of  these  endblocks 
were  separated  by  less  than  rc.  Next,  we  identified  the  con¬ 
formation  of  the  midblocks.  The  midblock  was  considered  to 
be  in  a  loop  conformation  if  both  endblocks  belong  to  the 
same  micelle,  whereas  if  the  midblock  connected  two  differ¬ 
ent  micelles,  then  it  was  considered  to  be  in  a  bridge  confor¬ 
mation.  Finally,  an  unassociated  endblock  also  was  assumed 
to  form  a  dangling  tail  conformation  for  the  midblock. 

To  explore  the  system  size  effects,  we  performed  equilibrium 
and  nonequilibrium  DPD  simulations  of  the  A1B6A1  triblock 
copolymer  in  midblock  selective  solvent  for  various  system 
sizes:  5000,  15,000,  and  30,000  DPD  particles,  where  results 
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FIGURE  4  (a)  Bridge  fraction,  and  (b)  equilibrium  modulus,  Ge  as  a  function  of  copolymer  concentration,  c/wt  %  from  the  exper¬ 
imental  measurements  of  the  SIS  system  (7)  and  DPD  simulations  of  the  A^A-i  triblock  copolymer  copolymer  for  simulated  sys¬ 
tem  sizes  of  5000,  15000,  and  30000  DPD  particles. 


are  presented  in  Figure  4.  From  Figure  4,  it  is  clear  that  the 
bridge  fraction  and  elastic  modulus  at  c  >  0.2  do  not  signifi¬ 
cantly  depend  on  the  system  size  (the  error  is  within  4%), 
but  the  bridge  fraction  and  elastic  modulus  at  c  —  0.2  are 
indeed  strongly  affected  by  the  system  size.  The  simulation 
results  for  15,000  and  30,000  particles  are  similar  for  the 
entire  concentration  range.  On  the  basis  of  these  findings, 
we  studied  the  morphology  and  rheology  of  the  triblock  co¬ 
polymer  in  the  midblock  selective  solvent  at  c  >  0.2  for  the 
entire  range  of  midblock  size  with  5000  total  number  of 
DPD  particles.  For  these  simulations,  we  used  a  periodic  box 
with  dimensions  11.9  x  11.9  x  11.9  (in  reduced  DPD  units), 
which  is  more  than  3x  larger  than  the  maximum  diameter 
of  the  micelles  obtained  from  our  simulation.  Therefore,  we 
do  not  expect  any  significant  finite  size  effects.  However,  we 
carried  out  equilibrium  and  nonequilibrium  simulations  at  c 
=  0.2  with  a  larger  system  of  15,000  DPD  particles  for  the 
entire  range  of  the  midblock  size.  In  addition,  we  have  car¬ 
ried  out  an  equilibrium  DPD  simulation  of  the  AxB^Ax  tri¬ 
block  copolymer  with  15000  DPD  particles  at  various  con¬ 
centrations  to  study  the  effect  of  concentration  and  relative 
block  size  on  the  morphology  of  the  system. 

RESULTS  AND  DISCUSSION 
Validation  of  Our  Method 

The  poly(styrene-isoprene-styrene)  block  copolymer,  or  SIS, 
in  an  isoprene-selective  solvent  has  been  chosen  as  a  model 
triblock  copolymer  for  this  study  because  of  the  sufficient 
amount  of  experimental  structural  and  rheological  data  avail¬ 
able.1'2'7'25  To  validate  our  method,  we  choose  A1B6A1  in  a  B- 
selective  solvent  as  a  prototype  of  the  SIS  copolymer  in  the 
I-selective  solvent,  n-tetradecane,  studied  by  Watanabe  et  al.7 
Molecular  weights  for  the  S-endblock  and  the  I-midblock  are 
7.2  kg/mol  and  35.8  kg/mol,  respectively.  Dynamics  mea¬ 
surements  have  shown  from  that  the  bead-spring  model 
accurately  describes  polystyrene  if  it  has  a  single  bead¬ 
spring  unit,  b  equal  to  5.0  nm,  corresponding  to  ~50  mono¬ 
mers.42'43  By  mapping  ~70  monomers  of  the  polystyrene 
block  to  a  single  DPD  particle,  we  anticipate  obtaining  a  rea¬ 


sonably  accurate  mesoscale  representation  of  our  system 
using  a  chain  of  8  DPD  particles. 

At  room  temperature,  the  Flory-Huggins  interaction  parame¬ 
ter  for  the  copolymer  blocks  Xps-pi  has  been  estimated  as 
0.1 144  (which  corresponds  to  yM  ~  90),  whereas  for  the 
polystyrene  block  and  n-tetradecane,  Xps-tetra  has  been  esti¬ 
mated  as  1.83  45  These  parameters  are  matched  to  our  co¬ 
polymer  model  in  a  strongly  selective  solvent  with  a  large 
degree  of  incompatibility  of  the  blocks. 

We  have  studied  the  following  weight-percent  concentrations 
of  AxBeAx,  0.2,  0.3,  0.4,  and  0.5,  which  correspond  to  volume 
fractions  0.17,  0.28,  0.39,  and  0.50,  respectively.  The  triblock 
is  expected  to  adopt  micellar  morphology  at  room  tempera¬ 
ture  for  this  concentration  range.  All  simulation  runs 
resulted  in  microphase  separation,  which  led  to  a  three- 
dimensional  network  structure  of  a  micellar  gel  with  the 
midblock  exhibiting  bridge  or  loop  conformations.  This  find¬ 
ing  is  in  agreement  with  the  experimental  study  of  Laurer 
et  al.1  For  a  similar  system,  SIS/mineral  oil,  they  observed 
micellar  morphology  at  a  styrene  concentration  <0.2  weight- 
percent.1  Figure  3  shows  an  example  of  an  equilibrated  sys¬ 
tem  from  our  DPD  simulations  for  c  =  0.2  and  c  =  0.4. 
Figure  3(a)  shows  a  separation  of  copolymer  solution  into 
densely  packed  micellar  phase  and  "soft"  (solvent)  phase  at 
c  =  0.2.  We  will  talk  about  this  phenomenon  in  details  in 
the  latter  sections. 

Figure  4(a)  compares  the  estimated  loop  fractions  of  the 
midblock  from  the  DPD  simulation  with  those  obtained  by 
Watanabe  et  al.7  by  measuring  the  dielectric  losses  of  the  SIS 
in  n-tetradecane.  Simulation  results  are  in  good  qualitative 
agreement  with  the  experimental  data  at  all  concentrations. 

The  storage  modulus,  G'  obtained  from  our  simulations  has  a 
numerical  value  close  to  a  plateau  for  the  frequencies  acces¬ 
sible  in  a  DPD  simulation.  (The  loss  modulus  G"[co )  is  not  a 
point  of  interest  in  this  study.)  The  equilibrium  value  of  the 
storage  modulus,  Ge  determined  from  our  DPD  simulations 
has  been  mapped  back  to  its  actual  value  using  eq  16.  Figure 
4(b)  shows  the  simulated  and  experimental  rheological  data 
and  illustrates  the  relatively  accuracy  of  the  DPD  model.  The 
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FIGURE  5  Bridge  fraction,  </>b  for  A-iB^A-i  copolymer  as  a  func¬ 
tion  of  the  relative  block  size,  x  at  various  volume  fractions  of 
the  triblock,  c. 

I-blocks  in  the  SIS  solutions  are  minimally  entangled  due  to 
their  low  molecular  weight  (Mw  «  35000],  when  compared 
with  entanglement  molecular  weight,  M'e  =  Me/c  with  Me  = 
6200  is  entanglement  molecular  weight  in  the  bulk.46  Hence, 
DPD  simulation  results  are  expected  to  be  in  good  agree¬ 
ment  with  the  experimental  data. 

Bridge  Fraction 

To  compare  the  effects  of  the  concentration  and  relative 
block  size  on  the  equilibrium  bridge  fraction,  <f>h,  we  have 
performed  a  set  of  equilibrium  DPD  simulations  with  vari¬ 
able  concentration,  c,  and  relative  volume  of  the  midblock,  x, 
where  results  are  presented  in  Figure  5.  Our  results  indicate 
that  0b  increases  with  increasing  copolymer  concentration; 
however,  increases  of  </>b  are  somewhat  less  significant  at 
higher  concentrations.  Likewise,  we  have  found  that  </>b 
increases  with  increasing  relative  volume  of  the  midblock  if 
the  value  of  the  midblock  volume  is  of  the  order  of  magni¬ 
tude  of  the  endblock  volume  (x  <  9].  For  the  relatively  large 
midblock  compared  with  the  endblock  (x  >  9],  </>b  only 
weakly  depends  on  x  and  reaches  some  constant  value  in  the 
range  of  0.40  to  0.66  depending  on  copolymer  concentration. 
This  behavior  is  a  result  of  the  reduction  of  the  entropic 
penalty  of  stretching  the  bridging  midblock.  A  decrease  in 
the  bridge  stretching  can  be  achieved  either  with  decreasing 
of  the  average  distance  between  the  micelles  as  copolymer 
concentration  increases,24  or  from  an  increase  in  the  relative 
size  of  the  midblock.  When  the  concentration  is  sufficiently 
high,  the  distance  between  micelles  becomes  shorter;  conse¬ 
quently,  the  effect  of  concentration  on  the  bridge  fraction  is 
weaker.  Decreasing  the  average  spacing  between  micelles 
with  increasing  concentration  has  also  been  reported  by 
Seitz  et  al 47  Likewise,  the  effect  of  the  relative  block  size  on 
the  bridge  fraction  becomes  negligible  when  the  midblock  is 
sufficiently  long,  so  that  micelles  bridge  without  significant 
stretching.  The  compact  size  of  the  midblock  compared  with 
the  distance  between  micelles  can  also  lead  to  the  tendency 
to  form  loops  for  low  x.5  For  example,  for  x  ==  3,  all  mid¬ 
blocks  form  loops  for  the  entire  concentration  range.  It  is 
important  to  mention  that  </>b  does  not  depend  on  the  abso¬ 
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lute  volume  of  the  blocks  but  rather  depends  on  the  relative 
block  sizes. 

For  a  relatively  high  degree  of  incompatibility  and  strongly 
selective  solvent  conditions,  0b  is  only  a  function  of  concen¬ 
tration  and  the  triblock  architecture.  However,  in  general,  0b 
also  depends  on  the  degree  of  incompatibility  and  miscibility 
of  the  copolymer  blocks  and  solvent,  0b  =  /  (c,  x,  yzy],  where 
indices  /  and  j  represent  A,  B,  and  S  particles.  This  depend¬ 
ence  becomes  significant  for  endblocks  with  a  low  degree  of 
polymerization  in  a  slightly  selective  solvent.  For  this  case, 
some  endblocks  would  not  associate  in  micelles,  producing 
dangling  tails  and  subsequently  decreasing  both  </>b  and  </>b24 
To  study  this  behavior,  we  could  systematically  vary  the  re¬ 
pulsive  parameters  in  our  model  and  then  evaluate  0b;  how¬ 
ever,  we  leave  such  a  study  for  future  work. 

Although  our  simulation  model  allows  reorganization  of  mid¬ 
blocks  during  the  simulation,  the  bridge-to-loop  ratio  can  be 
assumed  constant  for  the  length  of  our  simulation  runs.  To 
verify,  we  have  continued  to  run  the  A1B12A1  system  for 
another  six  simulation  cycles  (simulation  time  was  6000  in 
DPD  units]  and  found  that  the  average  bridge  fraction 
changed  less  than  1%.  Also,  we  have  checked  the  character¬ 
istic  time  of  the  bridge-to-loop  transformation  and  found 
that  it  is  a  relatively  long  event,  t  ~  103  in  DPD  units.  Figure 
6  shows  the  average  exchange  time  of  the  midblock  for  the 
AiB12A!  triblock  copolymer.  The  loops  become  less  stable  as 
concentration  increases.  This  finding  is  related  to  the 
entropic  penalty  associated  with  forcing  both  endblocks  to 
reside  in  the  same  micelle  47  Bridges  are  more  short-lived  at 
c  =  0.2  as  a  result  of  stretching  in  the  gel  exhibiting  synere- 
sis.  At  c  >  0.2,  the  bridge-to-loop  exchange  time  does  not 
depend  on  concentration. 

Syneresis  of  ABA  Gels 

According  to  our  simulation  results,  the  micelles  of  ABA 
copolymers  are  typically  formed  at  a  distance  comparable 
with  the  midblock  size  to  avoid  both  bridge  stretching.  If  the 
concentration  of  the  chains  is  small  and  the  midblock  is  rela¬ 
tively  short,  micelles  tend  to  be  crowded,  and  this 


FIGURE  6  The  characteristic  time  of  the  bridge-to-loop  and 
loop-to-bridge  exchange  as  a  function  of  concentration  for  the 
A-1B-12A-1  copolymer. 
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FIGURE  7  Vol  ume  fraction  of  the  solvent  phase,  f,  as  a  function  of  (a)  the  midblock  size,  x  at  c  =  0.2  and  (b)  copolymer 
3concentration,  c  for  the  A^A-i  and  A-iB^A-i  copolymers. 


occurrence  consequently  results  in  syneresis  wherein  separa¬ 
tion  of  the  polymer  solution  into  two  macrophases,  i.e.,  a 
densely  packed  micelle  phase  and  a  solvent-rich  phase.11'48 
Figure  3  [a]  shows  an  example  of  a  gel  exhibiting  syneresis 
with  a  distinct  “soft"  region  (solvent  phase]  in  a  copolymer 
solution  at  low  c.  We  have  probed  an  extent  of  syneresis  by 
measuring  the  volume  fraction  of  the  solvent  phase,  /  and 
defined  it  as  /  ==  where  Rfree  and  Vt  are  volume  of  sol¬ 
vent  phase  and  total  volume  of  the  system,  respectively.  We 
have  partitioned  our  simulation  box  into  bins,  where  the  bin 
size  is  approximately  equal  to  the  average  size  of  the 
micelles.  Then,  we  computed  the  number  of  the  bins  occu¬ 
pied  by  pure  solvent,  n0.  We  have  evaluated  /  from  /  = 
flol4>in,  where  Vbin  is  the  volume  of  the  bin.  We  have  pre¬ 
sented  our  results  in  Figure  7.  Our  simulations  have  shown 
that  /  increases  with  decreasing  concentration  and  block 
size.  For  the  copolymer  with  a  small  midblock,  the  syneresis 
occurs  even  for  higher  concentration,  for  example  for  the 
AiBbAx  triblock  at  c  —  0.3. 

On  the  basis  of  the  observation  of  the  system  undergoing 
syneresis,  we  expect  to  find  the  elastic  modulus  lower  at  c  = 
0.2  than  the  one  estimated  from  morphology  calculations. 
We  will  discuss  this  feature  in  the  latter  section. 


Elastic  Shear  Modulus 

To  assess  the  effect  of  concentration  and  relative  block  size 
on  the  equilibrium  elastic  modulus,  Ge,  we  have  performed 
oscillatory  shear  simulations  for  the  same  set  of  concentra¬ 
tions  and  relative  midblock  volumes  used  in  the  micelle 
bridging  study.  As  an  example  of  the  results  from  our  study, 
in  Figure  8,  we  present  oscillatory  shear  simulations  for  the 
AiBisAi  copolymer.  The  storage  modulus,  (7(cl>]  of  this  sys¬ 
tem  is  expected  to  obey  the  modified  Rouse  theory  for  net¬ 
worked  polymers,8  where  at  high  frequencies  G'[co ]  exhibits 
a  power  law  decay  with  exponent  0.53,  followed  by  a  plateau 
region,  corresponding  to  Ge  at  low  frequencies  (Fig.  8a].  The 
slope  of  the  G'[co)  decay  corresponds  to  a  bead-spring  model 
with  partial  hydrodynamic  interactions.  Thus,  we  have 
observed  the  intermediate  state  between  free-draining 
behavior  (where  chains  move  independently  of  each  other] 
and  hydrodynamic  behavior  (where  frictional  interactions 
dominate]  that  corresponds  to  1/2  and  2/3  power  laws, 
respectively, 8  Because  the  frequency  in  a  DPD  simulation  is 
expressed  as  \ where  MA  is  the  molecular  weight  of 
the  A-block,  the  decay  of  the  storage  modulus  is  nearly  inac¬ 
cessible  by  experiment  because  of  too  high  frequency  regime 
and,  therefore,  is  beyond  the  scope  of  our  study.  Because 
gels  are  characterized  by  G'}  which  are  nearly  independent  of 


FIGURE  8  Storage  modulus,  G  for  the  A-iB-ibA-i  triblock  copolymer  as  a  function  of  frequency,  co:  (a)  example  of  the  storage  modu¬ 
lus  for  concentration,  c  =  0.4  and  period  p  =  2II/co  from  10  to  10000.  The  solid  lines  are  approximations  using  a  quadratic  Bezier 
curve  representation,  (b)  The  plateau  region  of  the  storage  modulus,  G  at  various  concentrations.  Solid  lines  represent  linear 
extrapolations  for  accessing  experimental  frequencies. 
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FIGURE  9  Contribution  of  a  single  bridge,  Gb  and  loop,  G\  to 
the  equilibrium  elastic  modulus  as  a  function  of  the  relative 
block  size,  x  at  various  volume  fractions  of  the  copolymer,  c. 
Gb  and  G\  are  plotted  in  DPD  units.  The  arrow  indicates  G\  eval¬ 
uated  for  the  A-iBsA-i.  All  other  data  points  correspond  to  Gb. 

co  over  several  orders  of  magnitude,  we  can  exploit  the  pla¬ 
teau  region  of  G'  for  the  subsequent  analysis.  An  example  of 
the  plateau  regime  of  G'  as  a  function  of  concentration  is 
shown  in  Figure  8. 

Following  the  procedure  described  earlier  in  Viscoelastic 
Property  Calculations  Section,  we  have  obtained  Ge  for  vari¬ 
ous  concentrations  and  midblock  sizes  and  presented  our 
results  in  Figure  9.  From  our  bridge  fraction  study  in  Bridge 
Fraction  Section,  we  determined  that  the  midblock  of  A^Ai 
forms  only  loops;  therefore,  we  can  estimate  the  elastic  mod¬ 
ulus  of  a  single  loop  to  the  elastic  modulus,  G\  from  eq  15  by 
setting  </>b  equal  to  0.  Assuming  that  G\  does  not  change  with 
concentration  or  polymer  length,  we  have  estimated  G\  = 
0.11  (in  DPD  units).  By  using  this  quantity  in  eq  15,  we  have 
found  the  elastic  modulus  of  a  single  bridge,  Gb;  results  are 
shown  in  Figure  9.  Because  Gh  is  typically  an  order  of  magni¬ 
tude  larger  than  G\,  we  have  concluded  that  the  contribution 
of  the  loops  compared  with  the  bridges  is  negligible  for  the 
unentangled  polymeric  gel.  This  finding  contradicts  with  ex¬ 
perimental  observations  that  the  contributions  of  the  individ¬ 
ual  loops  and  tails  to  the  equilibrium  elastic  modulus  may 
become  comparable  with  bridges.7  The  explanation  of  this 
disagreement  is  the  following.  Because  of  the  soft  potentials 
we  have  used  in  this  study,  our  simulations  cannot  realisti¬ 
cally  describe  excluded  volume  effects.  As  a  result,  the  poly¬ 
mer  chains  can  pass  through  each  other  without  an  effective 
barrier  to  chain  crossing.  This  behavior  leads  to  a  diminished 
effect  on  the  loop  elasticity,  i.e.,  we  have  found  G\  to  be  an 
order  of  magnitude  smaller  than  the  measured  experimental 
value.7  For  the  simulation  sizes  used  in  this  study,  the 
domains  do  not  exhibit  long-range  order  in  their  spatial 
arrangement.  This  effect  could  also  contribute  to  the  reduc¬ 
tion  of  the  loop  elasticity.  Because  of  these  features,  we  can 
assume  that  the  loop  contribution,  Gh  does  not  depend  on 
concentration  nor  midblock  size  and  is  therefore  negligible. 
To  account  for  excluded  volume  and  entanglement  effects,  in 
future  work  we  will  improve  our  studies  by  introducing  a 


noncrossing  condition  within  the  DPD  method,49'50  which 
may  serve  to  mimic  the  entanglement  constraints  of  the  real 
systems. 

From  entropic  elasticity  theory,8  the  contribution  of  a  single 
bridge  to  the  equilibrium  elastic  modulus  can  be  estimated 
from 

Gb=gr\kBT,  (18) 

ro 

where  g  is  a  coefficient  on  the  order  of  unity,  determined  by 
a  method  for  calculating  chain  entropy.3  r\  is  the  mean 
square  end-to-end  distance  of  a  bridge,  and  r fj  is  the  mean 
square  end-to-end  distance  of  the  same  bridge,  if  it  is  not 
constrained  by  a  network.  (Note  that  kBT  is  1  for  our  simula¬ 
tion.)  r|  depends  on  the  average  distance  between  micelles 
and  characterizes  the  stretching  of  the  midblocks.47  Our  sim¬ 
ulations  have  found  that  the  average  distance  between 
micelles  decreases  with  concentration  (data  are  not  shown). 
The  ratio  r|/rb  would  be  equal  to  1  for  the  temperature  and 
concentration  at  which  a  network  was  produced  if  there  was 
no  midblock  stretching.8  However,  for  the  short  midblocks 
and  the  low  copolymer  concentration,  r\jr\  is  expected  to 
be  higher  than  unity.  Figure  9  reflects  stretching  of  the  mid¬ 
block  bridges  for  the  majority  of  our  data,  which  is  consist¬ 
ent  with  Gh  >  1.  At  c  —  0.2,  Gh  becomes  less  than  1  for  x  < 
15.  We  attribute  this  to  syneresis  that  has  occurred  (Fig.  3a), 
where  the  "soft"  solvent  phase  primarily  determines  the  elas¬ 
tic  modulus,  and  thus,  the  contribution  of  a  single  bridge  is 
diminished.  This  finding  is  in  agreement  with  experimental 
data  and  the  conclusions  of  Watanabe  et  al.7  Gh  also 
decreases  for  x  =  18  for  the  reason  that  these  long  mid¬ 
blocks  are  not  considerably  stretched.  This  effect  becomes 
more  prominent  for  c  =  0.5,  reflecting  the  average  decrease 
of  micelle  spacing  for  higher  copolymer  concentration.  We 
have  used  the  combined  variable  g  —  g  =  d(rE/ro )  by  aver¬ 
aging  our  data  for  all  concentration  and  block  sizes,  exclud¬ 
ing  Gb  >  1  ±  0.1  and  have  found  g  &  1.65.  This  parameter 
is  expected  to  characterize  approximately  the  equilibrium 
modulus  for  x  =  9  to  18. 

Mapping  our  results  to  actual  values  of  the  elastic  modulus 
for  the  ABA  triblock  copolymer,  we  have  finally  obtained  the 
following  conventional  equation  for  the  elastic  modulus  for 
unentangled  gels: 

Ge=d^c<Mx,c,Z(/)’  (19] 

In  eq  19,  the  volume  fraction,  c,  the  molar  volume  of  the  A- 
block,  vA,  and  the  degree  of  incompatibility,  Xij>  characterize 
the  physical  triblock  copolymer  and  solvent.  The  relative  size 
of  the  midblock  with  respect  to  the  endblock,  x,  is  evaluated 
from  eq  12.  Using  these  parameters,  the  bridge  fraction  0b 
(x,  c,  Xij )  is  determined  from  equilibrium  DPD  simulation. 
The  fit  from  this  analysis  of  the  equilibrium  shear  modulus 
and  the  simulation  results  are  shown  in  Figure  10,  where 
one  can  see  that  Ge  matches  well  with  simulation  results  for 
triblock  copolymers  at  c  >  0.2  with  large-sized  midblock 
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FIGURE  10  Elastic  modulus,  Ge  for  A-iB^A-i  as  a  function  of  the 
relative  block  size,  x  at  various  volume  fractions,  c.  Simulation 
data  are  represented  with  points  and  solid  lines  are  fit  to  eq  19 
with  [i  =  1 .65. 

relative  to  the  endblock  (x  >  9]  if  /i  is  equal  to  1.65.  Fur¬ 
thermore,  we  have  found  agreement  with  the  work  of 
Raspaud  et  al.  that  the  plateau  for  the  elastic  modulus  in  a 
selective  solvent  for  gels  in  the  semidilute  regime  made  from 
unentangled  chains  is  comparable  to  the  Rouse  modulus.9 

The  elastic  modulus  at  c  <  0.2  obtained  from  nonequilibrium 
simulation  is  significantly  lower  than  its  value  estimated  from 
the  bridge  fraction.  This  disagreement  between  results 
obtained  from  equilibrium  and  nonequilibrium  simulations  is 
explained  by  the  separation  of  liquid  from  the  gel  at  the  small 
concentration  [Fig.  3 (a}].  As  a  consequence,  we  cannot  predict 
the  modulus  at  c  =  0.2  from  equilibrium  simulation  alone. 

Our  simulations  show  the  elastic  modulus  of  unentangled  tri¬ 
block  chains  with  the  same  volume  of  the  endblock  for  each 
concentration  has  a  maximum  value  and  then  decreases  with 
increasing  x.  The  elastic  modulus  of  unentangled  polymeric 
gel  in  a  strongly  selective  solvent  seems  to  have  a  broad 
peak  for  x  between  6  and  12. 

Assuming  the  same  concentration  and  volume  of  the  triblock 
copolymer,  the  elastic  modulus  in  eq  19  depends  only  on  the 
bridge  fraction.  For  the  case  of  XabM  >  50,  </>b  increases  with 
increasing  relative  volume  of  the  midblock  for  x  <  9  and 
reaches  a  plateau  for  the  midblock  size  for  x  >  9,  which  cor¬ 
responds  to  /A  >  0.2  (eq  12]  for  each  copolymer  concentra¬ 
tion.  Taking  into  account  these  findings,  we  predict  that  the 
maximum  elastic  modulus  of  the  gel  made  of  ABA  triblock 
copolymers  of  low  molecular  weight  in  a  B-block  selective 
solvent  for  the  particular  polymer  size  and  concentration 
occurs  if  the  volume  fraction  of  the  endblock  in  the  "dry" 
copolymer  is  less  than  20%. 

The  main  advantage  of  this  approach  (eq  19]  is  that  it  can 
estimate  the  elastic  response  of  ABA  triblock  copolymers  of 
low  molecular  weight  in  a  midblock-selective  solvent  directly 
from  the  morphology  of  the  gels  obtained  from  equilibrium 
DPD  simulations  by  evaluating  (j)h  (x,  c,  Xij )■  Because  nonequi¬ 
librium  shear  simulation  is  not  necessary,  this  approach 
greatly  reduces  the  simulation  time. 


Equation  19  is  valid  only  for  block  copolymers  of  relatively 
low  molecular  weight  where  entanglement  does  not  likely 
occur.  For  copolymers  with  higher  molecular  weight,  entan¬ 
glements  would  give  rise  to  a  large  additional  contribution 
to  the  elastic  modulus.8 

CONCLUSIONS 

The  equilibrium  and  nonequilibrium  oscillatory-shear  DPD 
simulation  methods  have  been  used  for  mesoscale  simulation 
of  gel-forming  ABA  triblock  copolymer  in  a  midblock-selective 
solvent.  Our  results  indicate  that  the  bridge  fraction,  </>b, 
increases  with  increasing  copolymer  concentration  and  is  less 
significant  at  higher  concentrations.  We  have  demonstrated 
that  for  a  strongly  selective  solvent  and  for  XabM  >  50,  0b 
increases  with  increasing  relative  volume  of  the  midblock. 
When  the  size  of  the  midblock  becomes  significantly  larger 
compared  with  the  endblock  size  (x  >  9],  </>b  reaches  a  plateau 
in  the  range  of  0.40  to  0.66  as  the  copolymer  concentration 
increases.  We  have  observed  that  the  elastic  response  of  a  gel 
made  of  ABA  triblock  copolymers  of  low  molecular  weight  in  a 
B-block  selective  solvent  and  for  /ABM  >  50  is  comparable 
with  the  Rouse  modulus.  This  fact  allows  us  to  estimates  the 
absolute  value  of  the  elastic  modulus  of  unentangled  ABA  tri¬ 
block  copolymers  directly  from  the  morphology  of  its  micellar 
gel  by  using  equilibrium  DPD  simulations.  We  have 
demonstrated  that  our  simulation  results  are  in  good  qualita¬ 
tive  agreement  with  experimental  data  and  theoretical  predic¬ 
tions.  Moreover,  our  work  supports  the  use  of  DPD 
simulation  methods  to  study  mechanical  and  structural  prop¬ 
erties  of  various  polymeric  gels.  The  knowledge  derived  from 
the  types  of  mesoscale  simulations  performed  here  will  reduce 
the  time  and  effort  necessary  to  perform  experiments, 
which  will  eventually  allow  us  to  tune  the  mechanical  proper¬ 
ties  of  ABA  triblock  copolymer  gels  for  specific  applications. 
However,  to  extend  our  approach  and  predict  mechanical 
properties  of  triblock  copolymer  of  high  molecular 
weight,  entanglements  should  be  implemented  in  the  DPD 
simulation,49'50  or  other  simulation  techniques  should  be 
exploited,  such  as  the  "Sliplink”  models  adopted  for  block 
copolymers.51  We  defer  exploration  of  entangled  networks 
made  of  triblock  copolymers  of  higher  molecular  weight  to 
later  work. 
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